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ABSTRACT 

To  understand  the  radiative  impact  of  tropical  thin  cirrus  clouds,  the  frequency  of  occurrence  and  optical 
depths  of  these  clouds  have  been  derived.  “Thin”  cirrus  clouds  are  defined  here  as  being  those  that  are  not 
detected  by  the  operational  Moderate  Resolution  Imaging  Spectroradiometer  (MODIS)  cloud  mask,  cor¬ 
responding  to  an  optical  depth  value  of  approximately  0.3  or  smaller,  but  that  are  detectable  in  terms  of  the 
cirrus  reflectance  product  based  on  the  MODIS  1.375-yLtm  channel.  With  such  a  definition,  thin  cirrus  clouds 
were  present  in  more  than  40%  of  the  pixels  flagged  as  “clear  sky”  by  the  operational  MODIS  cloud  mask 
algorithm.  It  is  shown  that  these  thin  cirrus  clouds  are  frequently  observed  in  deep  convective  regions  in  the 
western  Pacific.  Thin  cirrus  optical  depths  were  derived  from  the  cirrus  reflectance  product.  Regions  of  sig¬ 
nificant  cloud  fraction  and  large  optical  depths  were  observed  in  the  Northern  Hemisphere  during  the  boreal 
spring  and  summer  and  moved  southward  during  the  boreal  autumn  and  winter.  The  radiative  effects  of 
tropical  thin  cirrus  clouds  were  studied  on  the  basis  of  the  retrieved  cirrus  optical  depths,  the  atmospheric 
profiles  derived  from  the  Atmospheric  Infrared  Sounder  (AIRS)  observations,  and  a  radiative  transfer  model 
in  conjunction  with  a  parameterization  of  ice  cloud  spectral  optical  properties.  To  understand  how  these 
clouds  regulate  the  radiation  field  in  the  atmosphere,  the  instantaneous  net  fluxes  at  the  top  of  the  atmosphere 
(TOA)  and  at  the  surface  were  calculated.  The  present  study  shows  positive  and  negative  net  forcings  at  the 
TO  A  and  at  the  surface,  respectively.  The  positive  (negative)  net  forcing  at  the  TOA  (surface)  is  due  to  the 
dominance  of  longwave  (shortwave)  forcing.  Both  the  TOA  and  surface  forcings  are  in  a  range  of  0-20  W  m~^, 
depending  on  the  optical  depths  of  thin  cirrus  clouds. 


1.  Introduction 

Numerous  studies  have  demonstrated  that  thin  cirrus 
clouds  are  frequently  present  near  the  tropical  tropo- 
pause  (Winker  and  Trepte  1998;  McFarquhar  et  al.  2000; 
Massie  et  al.  2002;  Pfister  et  al.  2001;  Dessler  et  al.  2006). 
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These  clouds  have  been  identified  as  one  of  the  sources 
of  uncertainty  in  the  study  of  the  earth’s  radiation 
budget  and  climate  (Liou  1986;  Lynch  et  al.  2002). 
Pfister  et  al.  (2001)  classified  cirrus  clouds  into  two  types 
in  terms  of  their  formation  mechanisms  based  on  the 
analysis  of  temperatures  along  back  trajectories  run  from 
these  clouds.  One  type  is  thin  quasi-laminar  wisps  cirrus 
associated  with  in  situ  formation  by  cooling  of  air  on 
a  synoptic  scale  or  mesoscale;  the  second  type  is  slightly 
lower,  thicker,  more  structured  cirrus  originating  from  re¬ 
cent  outflow  of  convection.  Using  the  Halogen  Occultation 
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Experiment  (HALOE)  observations,  Massie  et  al.  (2002) 
found  that  half  the  cirrus  clouds  observed  during  HALOE 
stemmed  from  convective  blowoff,  whereas  the  other 
half  were  associated  with  in  situ  formation. 

Because  of  the  high  frequency  of  occurrence  of  thin 
cirrus  clouds,  the  effect  of  these  clouds  on  the  earth’s 
radiation  budget  can  be  significant.  For  example,  these 
clouds,  located  high  in  the  atmosphere,  absorb  longwave 
radiation  but  emit  radiation  at  very  low  temperatures, 
producing  local  heating  by  a  few  degrees  per  day  (Jensen 
et  al.  1996;  McFarquhar  et  al.  2000)  and  net  positive 
cloud  radiative  forcing  on  the  order  of  1  W  (Wang 
et  al.  1996;  McFarquhar  et  al.  2000).  Because  of  the 
difficulties  associated  with  in  situ  measurements,  it  is 
challenging  to  accurately  quantify  the  radiative  forcing 
of  these  clouds.  To  better  understand  the  effect  of  thin 
cirrus  clouds  on  climate,  there  is  a  pressing  need  to  im¬ 
prove  the  current  knowledge  of  the  radiative  effect  of 
thin  cirrus  clouds  near  the  tropopause. 

Several  studies  based  on  lidar  measurements  (Winker 
and  Trepte  1998;  Pfister  et  al.  2001;  Comstock  et  al. 
2002)  and  ground-based  measurements  (Comstock  and 
Sassen  2001;  Mather  et  al.  1998;  Platt  et  al.  1998)  pro¬ 
vided  valuable  results,  but  with  temporal  and  spatial 
sampling  limitations.  The  intent  of  this  study  is  to  un¬ 
derstand  the  climatological  characteristics  of  tropical 
thin  cirrus  optical  depths,  in  particular  those  that  are  not 
detected  by  imagers  using  standard  methods,  and  to 
quantify  the  radiative  forcing  of  these  clouds. 

2.  Data 

A  calculation  of  cloud  radiative  forcing  requires 
cloud  optical  properties  as  well  as  the  profiles  of  at¬ 
mospheric  water  vapor  and  temperature.  In  this  study, 
thin  cirrus  optical  properties  are  derived  from  the  op¬ 
erational  Moderate  Resolution  Imaging  Spectroradi- 
ometer  (MODIS)  cirrus  reflectance  dataset  from  the 
Earth  Observing  System  (EOS)  Aqua  satellite.  MODIS 
has  36  distinct  spectral  bands  (20  reflected  solar  bands 
and  16  thermal  emissive  bands)  with  wavelengths  rang¬ 
ing  from  0.41  to  14.2  jum.  In  this  study,  level-2  cirrus 
reflectances  (contained  in  the  Aqua  MODIS  product 
MYD06)  derived  from  the  MODIS  0.66-  and  1.375-/im 
band  measurements  (Gao  et  al.  2002)  during  a  period 
from  June  2005  to  May  2006  are  used  to  retrieve  the 
corresponding  cirrus  optical  depths.  The  retrieval  is  per¬ 
formed  at  every  pixel  between  30°S  and  30°N  latitudes, 
with  a  1-km  horizontal  resolution.  Additionally,  level-2 
cloud  mask  data  (MYD35)  are  used  to  limit  cirrus  op¬ 
tical  depth  retrievals  to  those  pixels  identified  as  “cloud 
free”  by  the  cloud  mask.  Recent  studies  suggest  that  the 
MODIS  cloud  mask  has  limited  detection  capability  when 


cirrus  optical  depths  are  less  than  about  0.3  (Ackerman 
et  al.  2008). 

Measurements  of  tropospheric  specific  humidity  (q) 
and  temperature  (T)  are  obtained  from  measurements 
made  by  the  Atmospheric  Infrared  Sounder  (AIRS)  on 
Aqua  (Aumann  et  al.  2003).  The  q  and  T  data  with  an 
almost  global  coverage  (Fetzer  et  al.  2006)  have  vertical 
resolutions  of  2  and  1  km,  respectively,  and  accuracies 
of  10%  and  1  K,  respectively.  Sampling  biases  in  the 
AIRS  datasets  are  minimized  by  the  retrieval  ability  of 
the  instrument  in  the  presence  of  up  to  80%  cloud  cover 
(Susskind  et  al.  2003).  The  horizontal  resolution  of 
a  single  retrieval  is  40-50  km. 


3.  Methodology 

a.  Optical  depth  of  tropical  thin  cirrus  clouds 

Cirrus  clouds  are  often  present  in  the  upper  tropo¬ 
sphere  or  lower  stratosphere,  and  more  than  90%  of 
the  total  water  vapor  in  the  atmosphere  is  distributed 
beneath  them  (Gao  and  Kaufman  1995).  In  relatively 
humid  environments,  the  MODIS  1.375-/im  channel  re¬ 
ceives  very  little  solar  radiance  reflected  by  the  surface 
or  the  lower  portion  of  the  atmosphere  (including  low- 
level  water  clouds)  when  high-level  cirrus  clouds  are  not 
present.  This  is  because  the  radiation  within  this  band  is 
absorbed  by  water  vapor  in  the  lower  atmosphere  under 
cirrus-free  conditions.  However,  if  high  clouds  are  pres¬ 
ent,  strong  signals  associated  with  the  reflection  and 
scattering  of  the  incident  solar  radiation  by  these  clouds 
are  observed  by  the  MODIS  sensor  at  the  131 5- pm 
channel.  Note  that  the  reflected  solar  radiation  observed 
by  the  MODIS  sensor  is  slightly  attenuated  by  water 
vapor  above  the  clouds.  The  operational  MODIS  cirrus 
reflectance  data  product  adjusts  for  this,  determining  the 
reflectance  that  would  exist  in  the  absence  of  water 
vapor  absorption  (Gao  et  al.  2002). 

This  study  focuses  on  tropical  thin  and  subvisual  cirrus 
clouds  whose  optical  depths  are  quite  small  (<Cl).  In  this 
case,  the  first-order  scattering  approximation  is  valid  and 
the  optical  depth  (r)  can  be  derived  from  the  cirrus  re¬ 
flectance  rc  as  follows  (Dessler  and  Yang  2003): 


T  =  r  - - - 


(1) 


where  p  is  the  cosine  of  the  sensor  zenith  angle,  po  is  the 
cosine  of  the  solar  zenith  angle,  4>  is  the  sensor  azimuthal 
angle,  <po  is  the  solar  azimuthal  angle,  w  is  the  single¬ 
scattering  albedo,  and  P  is  the  scattering  phase  function. 
The  uncertainties  of  the  retrieved  optical  thickness 
based  on  Eq.  (1)  due  to  the  uncertainties  in  the  phase 
function  have  been  estimated  by  Dessler  and  Yang  (2003). 
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Fig.  1.  Fraction  of  clear-sky  observations  for  1°  X  1°  boxes  that  have  detectible  thin  cirrus  optical  depths  larger  than 
0.02  for  boreal  (a)  spring,  (b)  summer,  (c)  autumn,  and  (d)  winter.  The  definition  of  season  is  determined  by  con¬ 
ventional  boreal  season  such  that  spring  months  are  March,  April,  and  May;  summer  months  are  June,  July,  and 
August;  autumn  months  are  September,  October,  and  November;  and  winter  months  are  December,  January,  and 
February. 


A  method  similar  to  the  present  approach  has  been 
developed  by  Meyer  et  al.  (2004)  using  a  rigorous  radi¬ 
ative  transfer  model  to  generate  lookup  libraries  for 
the  T-Vc  relation,  which  includes  the  effect  of  multiple 
scattering  and  thus  is  applicable  to  both  thin  and  thick 
cirrus  clouds.  However,  the  present  method  is  more 
straightforward  and  computationally  efficient  although 
it  is  limited  to  thin  and  subvisual  cirrus  clouds.  It  should 
be  noted  that  various  complementary  techniques  have 


also  been  developed  in  the  past  two  decades  to  derive 
cirrus  microphysical  and  optical  properties  (e.g.,  Minnis 
et  al.  1993a,b,  1995;  King  et  al.  2004;  Kokhanovsky  and 
Nauss  2005;  Chepfer  et  al.  1998, 2002;  Stubenrauch  et  al. 
1999),  including  the  operational  MODIS  optical  and 
microphysical  cloud  product  (Platnick  et  al.  2003). 

Ice  crystals  in  cirrus  clouds  are  nonspherical  particles 
with  complex  shapes  (e.g.,  Heymsfield  1975;  Heymsfield 
and  laquinta  2000).  In  the  present  retrieval,  we  assume 
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that  a  cirrus  cloud  consists  of  a  mixture  of  bullet  rosettes, 
hollow  columns,  aggregates,  and  hexagonal  plates  with 
percentages  as  follows:  25%  plates,  25%  hollow  col¬ 
umns,  and  50%  bullet  rosettes  for  size  bins  if  the  particle 
maximum  dimension  is  smaller  than  70  /im;  and  20% 
plates,  20%  hollow  columns,  30%  bullet  rosettes,  and 
30%  aggregates  for  size  bins  larger  than  70  /mm  (Baum 
et  al.  2000;  Platnick  et  al.  2003;  King  et  al.  2004).  Spe¬ 
cifically,  the  phase  function  used  for  the  present  retrieval 
of  cirrus  optical  depth  corresponds  to  an  effective  par¬ 
ticle  radius  of  6.7  jum  [see  Fig.  4  in  King  et  al.  (2004); 
also  see  Karcher  and  Lohmann  (2002)  and  Garrett  et  al. 
(2003)  for  the  size  distribution  of  cirrus  clouds],  typical 
for  cold  cirrus  clouds.  The  average  cloud  top  pressure 
and  temperature  for  the  tropical  thin  cirrus  clouds  are 
200  mb  and  200  K  (Pfister  et  al.  2001). 

As  noticed  by  Dessler  and  Yang  (2003),  the  signals 
associated  with  the  reflection  of  thin  cirrus  are  quite 
weak.  Thus,  the  reflection  from  a  bright  underlying  sur¬ 
face  under  dry  atmospheric  conditions,  such  as  over  des¬ 
erts  where  the  radiation  reflected  by  the  surface  is  not 
completely  absorbed  by  water  vapor,  may  contaminate 
the  signals  associated  with  true  cirrus  reflection.  For  this 
reason,  this  study  is  confined  to  the  cloud-free  pixels 
over  the  oceans  where  the  surface  reflectance  is  rela¬ 
tively  homogeneous  and  small.  Dessler  and  Yang  (2003) 
also  analyzed  various  error  sources  and  determined  that 
the  detection  limit  in  terms  of  optical  depth  by  this 
method  is  0.02.  If  the  retrieved  optical  depth  is  larger 
than  0.02  for  a  pixel,  a  thin  cirrus  cloud  is  assumed  to  be 
present  within  this  pixel. 

b.  Radiative  transfer  model 

To  investigate  the  radiative  forcing  of  thin  cirrus 
clouds,  the  shortwave  and  longwave  fluxes  at  the  top  of 
the  atmosphere  and  at  the  surface  are  calculated  from 
a  radiative  transfer  model,  namely  LibRadtran  de¬ 
veloped  by  Mayer  and  Kylling  (2005).  The  absorption  of 
H2O,  O2, 03,  CO2,  and  N2O  is  fully  taken  into  account  in 
the  LibRadtran  code.  The  profiles  of  temperature,  H2O, 
and  O3  are  obtained  from  the  AIRS  level-3  monthly 
products.  Using  these  profiles,  background  gaseous  ab¬ 
sorption  properties  are  calculated  using  the  correlated 
k-distribution  method  (Fu  and  Liou  1992,  and  references 
therein).  Because  the  AIRS  products  do  not  provide  at¬ 
mospheric  profiles  at  altitudes  corresponding  to  pressures 
lower  than  100  hPa  for  water  vapor  and  1  hPa  for  tem¬ 
perature,  standard  tropical  atmospheric  profiles  are  used 
above  the  highest  AIRS  levels.  The  distributions  of  other 
gases  such  as  O2,  CO2,  and  NO2  are  obtained  from  the 
standard  tropical  atmospheric  profile.  These  profiles  are 
interpolated  into  100  layers  in  the  model  and  cloud  layer 
is  inserted  as  one  layer  from  14  to  15  km. 
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Fig.  2.  (a)  Histograms  of  optical  depth  of  thin  cirrus  clouds  re¬ 
trieved  between  latitudes  30°S  and  30°N  for  each  of  the  boreal 
seasons,  (b)  Cumulative  fractions  of  optical  depth  of  thin  cirrus 
clouds  for  each  of  the  boreal  seasons. 

c.  Bulk  optical  properties  of  ice  clouds 

The  radiative  properties  of  cirrus  clouds  are  funda¬ 
mentally  determined  by  their  microphysical  and  optical 
properties  such  as  ice  crystal  shape,  size  distribution,  ice 
water  content  (IWC),  and  effective  particle  radius  (R^). 
The  single-scattering  properties  of  ice  clouds,  namely 
the  asymmetry  parameter,  extinction  coefficient,  and 
single-scattering  albedo,  are  usually  parameterized  in 
terms  of  IWC  and 

In  the  past  three  decades,  significant  progress  has  been 
made  in  the  parameterization  of  the  bulk  optical  proper¬ 
ties  of  ice  clouds  (e.g.,  Ebert  and  Curry  1992;  Fu  1996;  Fu 
et  al.  1998;  Key  et  al.  2002;  McFarquhar  et  al.  2002; 
Edwards  et  al.  2007).  The  parameterization  of  cirrus  op¬ 
tical  properties  in  this  study  is  based  on  that  reported  by 
Key  et  al.  (2002)  but  with  an  extension  to  the  infrared 
spectral  band.  Similar  to  Key  et  al.  (2002),  thirty  size  dis¬ 
tributions  reported  by  Fu  (1996)  and  Mitchell  et  al.  (1996) 
are  used.  The  single-scattering  properties  of  individual  ice 
crystals  are  taken  from  Yang  et  al.  (2000,  2005). 
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Fig.  3.  Optical  depth  of  tropical  thin  cirrus  clouds  for  the  pixels  flagged  as  clear  sky  by  MODIS  for  boreal  (a)  spring, 

(b)  summer,  (c)  autumn,  and  (d)  winter. 


For  computational  efficiency  in  the  simulation  of  ra¬ 
diation  flux,  a  common  approach  is  to  approximate  the 
scattering  phase  function  in  terms  of  the  Henyey- 
Greenstein  function  (Henyey  and  Greenstein  1941). 
However,  the  Henyey-Greenstein  function  does  not  rep¬ 
resent  a  good  approximation  to  the  true  phase  function 
in  backscattering  directions,  particularly  in  the  case  of 
the  phase  functions  of  ice  crystals.  To  overcome  this 
shortcoming,  in  this  study  we  follow  Key  et  al.  (2002) 
and  use  the  double  Henyey-Greenstein  function.  A 
method  for  determining  the  parameters  involved  in  the 


double  Henyey-Greenstein  function  has  been  reported 
by  Key  et  al.  (2002). 

4.  Results 

a.  The  temporal  and  spatial  distribution  of  tropical 
thin  cirrus  clouds 

The  temporal  and  spatial  distributions  of  the  mean 
optical  depth  of  thin  cirrus  clouds  and  the  frequency  of 
occurrence  are  analyzed  for  a  period  from  June  2005 
to  May  2006.  The  analysis  is  designed  to  determine  the 
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locations  of  the  maximum  and  minimum  occurrences  of 
thin  cirrus  clouds  and  to  examine  the  characteristics  of 
their  seasonal  variations. 

Figure  1  shows  the  fraction  of  the  MODIS  observa¬ 
tions  made  for  clear-sky  conditions  according  to  the 
MODIS  cloud  mask,  when  optical  depths  obtained  from 
Eq.  (1)  are  larger  than  0.02.  Using  ground-based  lidar 
observation,  Comstock  et  al.  (2002)  showed  that  the 
median  optical  depth  of  single  layer  cirrus  clouds  whose 
cloud-base  heights  are  higher  than  15  km  is  0.019;  that  is, 
50%  of  the  optical  depths  of  single-layer  cirrus  clouds 
measured  by  ground-based  lidar  are  lower  than  0.019.  It 
is  evident  from  Fig.  1  that  high  fraction  occurs  in  the 
Northern  Hemisphere  during  the  boreal  spring  and 
summer  and  in  the  Southern  Hemisphere  during  the 
boreal  autumn  and  winter.  In  spring,  the  high  fraction 
occurs  over  the  western  Pacific  around  Indonesia.  The 
fraction  is  larger  than  50%  over  most  of  the  Northern 
Hemisphere  oceans  whereas  the  fraction  is  smaller  than 
30%  over  the  Southern  Hemisphere  oceans.  Over  the 
western  Pacific,  the  high  fraction  is  observed  around 
deep  convection  regions  where  thin  cirrus  clouds  may  be 
related  to  outflow  from  convective  anvil  cirrus  or  to  in 
situ  formation  in  the  cold  phase  of  waves  spawned  by  the 
convection. 

The  fraction  of  observations  for  summer  shows  spatial 
distributions  that  are  similar  to  those  for  spring.  High 
fraction  is  observed  over  the  western  Pacific,  the  Persian 
Gulf,  and  Central  America;  low  fraction  occurs  over  the 
southeastern  Pacific  and  the  southern  Indian  Ocean. 
High  fraction  occurrences  are  still  observed  over  the 
western  Pacific  but  move  southward  in  autumn.  The 
Southern  Hemisphere  oceans  show  high  fraction  whereas 
the  Northern  Hemisphere  oceans  show  less  than  30% 
occurrence.  Other  high  fraction  regions  are  off  the  west 
coast  of  South  America  and  off  the  west  coast  of  southern 
Africa.  The  high  fraction  occurs  over  the  western  Pacific 
during  winter  and  autumn  and  the  low  fraction  is  widely 
observed  over  Northern  Hemisphere  oceans. 

Figure  2a  shows  the  histograms  of  cirrus  optical 
depths  for  each  season.  The  seasonal  fraction  is  calcu¬ 
lated  with  a  resolution  of  0.005  for  cirrus  optical  depth. 
The  distributions  for  the  four  seasons  are  similar.  For 
a  given  season,  the  fraction  increases  sharply  to  its 
maximum  at  optical  depths  near  0.01,  decreases  rapidly 
to  small  values  at  optical  depths  near  0.04,  and  further 
decreases  at  large  optical  depths  to  quite  small  values. 
Figure  2b  shows  the  cumulative  fraction  of  optical  depth 
for  each  season.  Using  r  =  0.02  as  a  detection  limit  as 
explained  in  section  3,  44%,  39%,  47%,  and  44%  of  the 
observations  flagged  as  cloud-free  by  the  MODIS  op¬ 
erational  cloud  mask  algorithm  have  detectible  thin 
cirrus  clouds  during  the  boreal  spring,  summer,  autumn. 


Fig.  4.  Zonally  averaged  (a)  fraction  of  observations  and 
(b)  optical  depth  of  thin  cirrus  clouds  for  each  season. 


and  winter,  respectively.  Dessler  and  Yang  (2003)  showed 
that  the  percentage  of  detectible  thin  cirrus  was  28% 
and  38%  during  the  6-8  June  2001  and  6-8  December 
2000  periods,  respectively.  The  percentage  in  the  current 
study  increased  compared  to  that  of  Dessler  and  Yang 
(2003)  during  both  seasons.  The  difference  in  the  per¬ 
centage  between  the  two  studies  might  result  from  the 
difference  in  the  amount  of  analyzed  data  or  from  the 
fact  that  cloud  formation  increases  in  the  afternoon. 
Note  that  the  Terra  MODIS  data  used  in  the  study  by 
Dessler  and  Yang  (2003)  were  associated  with  morning 
properties  whereas  the  Aqua  MODIS  data  reflect  af¬ 
ternoon  properties.  Using  data  from  Terra  and  Aqua  for 
the  same  period  might  be  useful  to  understand  the  dif¬ 
ference  between  morning  and  afternoon  characteristics 
of  thin  cirrus  clouds. 

Figure  3  shows  the  spatial  distributions  of  the  average 
optical  depths  of  tropical  thin  cirrus  clouds  for  the  pixels 
flagged  as  clear  sky  for  each  of  the  four  seasons.  The 
pattern  is  very  similar  to  that  of  the  observed  fraction. 
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Fig.  5.  Spatial  distributions  of  (a)  shortwave,  (b)  longwave,  and  (c)  net  cloud  radiative  forcing  at  the  top  of  the 

atmosphere  for  June  2005.  Units  are  W  m“^. 


During  spring,  large  optical  depths  are  observed  over 
the  western  and  northern  Pacific  and  small  optical 
depths  are  observed  over  the  Southern  Hemisphere 
oceans.  In  summer,  large  optical  depths  appear  widely 
over  the  western  Pacific,  the  northern  Indian  Ocean,  and 
Central  America  whereas  small  optical  depths  are  ob¬ 
served  over  the  Southern  Hemisphere  ocean  areas.  The 
locations  of  large  optical  depths  are  still  observed  over 
the  western  Pacific  but  migrate  southward  during  au¬ 
tumn.  Large  optical  depths  also  appear  at  the  west  coast 
of  South  America.  The  pattern  of  spatial  distributions  of 
optical  depth  for  winter  is  similar  to  that  for  autumn 
when  large  optical  depths  are  observed  over  the  western 
Pacific  and  off  the  west  coast  of  Africa. 

Figure  4  shows  the  zonally  averaged  fraction  of  ob¬ 
servations  and  optical  depths.  The  fraction  of  observa¬ 
tions  and  optical  depths  are  averaged  with  a  resolution 
of  1°  over  latitudes  from  30°S  to  30°N.  From  the  figure, 
latitudinal  variations  of  the  fraction  of  observations  and 
optical  depths  with  the  seasons  are  evident.  The  high 


fraction  and  large  optical  depths  are  observed  in  the 
Northern  Hemisphere  high  latitudes  and  decrease  to¬ 
ward  the  south  during  spring  and  summer.  On  the  other 
hand,  the  high  fraction  and  large  optical  depth  appear  in 
the  Southern  Hemisphere  at  high  latitudes  and  decrease 
toward  the  north  during  autumn  and  winter. 

h.  The  radiative  forcing  of  thin  cirrus  clouds 

Cloud  radiative  forcing  has  been  investigated  on  the 
basis  of  satellite  or  ground-based  measurements  (e.g., 
Ramanathan  et  al.  1989;  Rossow  and  Zhang  1995;  Mace 
et  al.  2006).  Clouds  have  a  warming  effect  on  the  surface 
or  in  the  atmosphere  by  absorption  and  emission  of 
longwave  radiation  and  a  cooling  effect  by  reflecting  the 
incident  shortwave  radiation.  The  net  cloud  radiative 
forcing  is  determined  by  the  sum  of  longwave  warming 
effect  and  shortwave  cooling  effect.  The  cloud  radiative 
forcing  is  defined  as  the  difference  between  clear  con¬ 
ditions  and  all-sky  conditions  (Ramanathan  et  al.  1989; 
Hartmann  et  al.  2001;  Yang  et  al.  2007). 
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Fig.  6.  As  in  Fig.  5,  but  at  the  surface  for  June  2005. 


The  solar  zenith  angle  from  each  pixel  is  averaged 
over  1°  X  1°  latitude  and  longitude  boxes  for  every 
month  and  is  incorporated  into  the  model  calculation. 
Clear-sky  fraction  of  observations  with  optical  depths 
larger  than  0.02  is  used  in  the  calculation  of  all-sky  flux. 
Flux  is  calculated  at  each  1°  X  1°  box  for  each  month 
using  the  monthly  averaged  optical  depth,  clear-sky 
fraction,  solar  zenith  angle,  and  corresponding  AIRS 
atmospheric  profile. 

Figure  5  shows  the  1°  X  1°  spatial  distribution  of  in¬ 
stantaneous  cloud  radiative  forcing  at  the  top  of  the  at¬ 
mosphere  calculated  from  the  retrieved  optical  depths, 
AIRS  profile,  and  averaged  solar  zenith  angle  for  June 
2005.  The  large  cloud  radiative  forcing  appears  around 
the  western  Pacific  where  the  optical  depths  are  large.  It 
shows  a  negative  cloud  radiative  forcing  for  shortwave 
radiation,  a  positive  forcing  for  longwave  radiation,  and 
a  net  positive  cloud  radiative  forcing. 

Figure  6  shows  the  spatial  distribution  of  instantaneous 
cloud  radiative  forcing  at  the  surface  for  June  2005.  Thin 
cirrus  clouds  produce  large  negative  shortwave  forcing  in 


the  Indian  Ocean  and  the  western  Pacific  Ocean  but  very 
small  positive  longwave  forcing  at  the  surface.  Thus,  the 
net  cloud  forcing  shows  a  negative  effect. 

Figure  7  shows  the  relationship  between  optical  depth 
and  instantaneous  cloud  radiative  forcing  with  respect  to 
the  solar  zenith  angle  obtained  from  1°  X  1°  calculations 
at  the  top  of  the  atmosphere  and  at  the  surface.  At  the 
top  of  the  atmosphere,  both  shortwave  negative  forcing 
and  longwave  positive  forcing  increase  in  magnitude 
with  increasing  optical  depths.  The  net  cloud  radiative 
forcing  shows  positive  values  of  0-10  W  because  of 
larger  positive  longwave  forcing.  The  difference  in  sur¬ 
face  temperature  and  atmospheric  profile  causes  the 
variation  in  the  radiative  forcing  when  the  optical  depth 
is  constant. 

At  the  surface,  the  radiative  forcing  of  thin  cirrus 
clouds  is  negative  and  positive  for  shortwave  and  long¬ 
wave  radiation,  respectively,  and  the  corresponding  net 
forcing  is  in  a  range  of  —20  to  0  W  m“^.  Shortwave 
negative  forcing  dominates  at  the  surface.  For  overcast 
sky  when  the  cloud  fraction  (N)  is  1  (not  shown  here). 
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Fig.  7.  (top)  Shortwave,  (middle)  longwave,  and  (bottom)  net  cloud  radiative  forcing  as  a  function  of  optical  depth  of 
thin  cirrus  clouds  and  solar  zenith  angle  (a)  at  the  top  of  the  atmosphere  and  (b)  at  the  surface. 


the  shortwave  negative  forcing  increases  in  magnitude 
with  the  increasing  solar  zenith  angle  if  the  optical  depth 
is  the  same.  For  longwave,  the  forcing  does  not  depend 
on  solar  zenith  angle.  However,  such  a  feature  is  not 
obvious  after  cloud  fraction  is  incorporated  into  the 
radiative  forcing  calculation. 

Table  1  lists  the  values  of  the  seasonally  averaged 
cloud  radiative  forcing  between  30°S  and  30°N.  Thin 
cirrus  clouds  induce  about  1  W  m“^  net  warming  at  the 
top  of  the  atmosphere  and  1.35  W  net  cooling  at 
the  surface.  The  possible  contamination  of  clear-sky 
pixels  by  thin  cirrus  clouds  may  result  in  inaccuracies  in 
the  estimation  of  fluxes  and  atmospheric  property  re¬ 
trievals.  Further  investigation  is  necessary  to  better  un¬ 
derstand  those  effects. 

5.  Summary 

In  this  study,  the  frequency  of  occurrence  of  tropical 
thin  cirrus  clouds  (defined  here  as  those  not  detected  by 
the  operational  MODIS  cloud  mask,  corresponding  to 
an  approximate  optical  depth  of  about  0.3  or  smaller), 
the  mean  optical  depth,  and  the  radiative  impact  were 
investigated.  The  optical  depth  of  thin  cirrus  clouds  were 
retrieved  using  the  MODIS  level-2  cirrus  reflectance 
product.  The  retrieval  was  confined  to  clear-sky  pixels 
over  the  oceans  to  remove  the  effect  from  the  surface  or 
lower  atmosphere.  Using  r  =  0.02  as  a  detection  limit, 
over  40%  of  the  observations  flagged  as  cloud-free  by 
MODIS  have  detectible  thin  cirrus  clouds.  The  regions 
of  high  occurrence  and  large  optical  depth  are  located 
around  deep  convection  areas  and  show  seasonal  vari¬ 


ations.  The  thin  cirrus  clouds  occur  more  frequently  with 
larger  optical  depths  in  the  Northern  Hemisphere  dur¬ 
ing  spring  and  summer  whereas  they  occur  more  fre¬ 
quently  in  the  Southern  Hemisphere  during  autumn  and 
winter. 

The  instantaneous  cloud  radiative  forcings  at  the  top 
of  atmosphere  and  at  the  surface  are  calculated  to  ex¬ 
amine  the  radiative  effect  of  this  subset  of  thin  cirrus 
clouds.  The  retrieved  optical  depth,  fraction  of  obser¬ 
vations,  atmospheric  profiles  from  AIRS,  solar  zenith 
angle,  and  parameterized  single-scattering  properties  of 
ice  clouds  are  used  to  simulate  the  radiative  flux  from  the 
LibRadtran  radiative  transfer  code.  Thin  cirrus  clouds 
have  a  net  positive  radiative  forcing  at  the  top  of  atmo¬ 
sphere  because  of  a  larger  positive  longwave  warming 
and  a  net  negative  radiative  forcing  at  the  surface  be¬ 
cause  of  a  larger  negative  shortwave  cooling.  To  improve 
current  understanding  of  the  properties  of  tropical  thin 
cirrus  clouds,  accurate  information  about  ice  crystal 
particle  size  and  cloud  height  may  be  necessary. 


Table  1.  The  seasonal  average  of  cloud  radiative  forcing  be¬ 
tween  30°S  and  30°N  latitudes.  All  units  are  W  ni“^.  (TOA  is  top  of 
atmosphere;  BOA  is  bottom  of  atmosphere.) 


Spring 

Summer 

Autumn 

Winter 

TOA  SW 

-1.74 

-1.62 

-1.79 

-1.69 

TOALW 

2.81 

2.56 

2.79 

2.69 

TOA  NET 

1.07 

0.94 

0.99 

1.00 

BOA  SW 

-2.01 

-1.86 

-2.06 

-1.94 

BOALW 

0.64 

0.59 

0.64 

0.60 

BOA  NET 

-1.37 

-1.28 

-1.41 

-1.34 
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